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Abstract 

For a class of tight-binding many-electron models on hyper-cubic lattices the 
equal-time correlation functions at non-zero temperature are proved to decay 
exponentially in the distance between the center of positions of the electrons 
and the center of positions of the holes. The decay bounds hold in any space 
dimension in the thermodynamic limit if the interaction is sufficiently small 
depending on temperature. The proof is based on the (7(l)-invariance prop- 
erty and volume- independent perturbative bounds of the finite dimensional 
Grassmann integrals formulating the correlation functions. 

1 Introduction 

A Hamiltonian governing the total energy of many electrons hopping and interact- 
ing on a finite lattice is defined as a self-adjoint operator on the finite dimensional 
Hilbert space of all states of electrons. Correlation functions in the system at non- 
zero temperature are formulated as a quotient of trace operations over the Hilbert 
space. Despite the explicitness of their mathematical definitions, to rigorously an- 
alyze the behavior of the correlation functions still requires restrictive assumptions 
and remains to be solved in a general setting. 

The method using Bogoliubov's inequality initiated by Hohenberg |^, Mermin 
and Wagner [12] has been applied to prove decay of various order parameters in the 
one- and two-dimensional Hubbard models. It has been shown that the magnetic 
order parameters ( [H] , [S] , [20] ) and the electron-pairing order parameters ([IH],[in]) 
vanish in the thermodynamic limit as the amplitude of the corresponding exter- 
nal field goes to zero. In these theories the application of Bogoliubov's inequality 
demonstrates that the thermodynamic limit of the order parameter under consid- 
eration is bounded from above by a constant times the inverse of an integral of the 
form 
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with the space dimension d and the amplitude A > of the external field. If this 
integral diverges as the amplitude A is sent to zero, which is true for d — 1,2, 
not true for d > 3, one concludes that the order parameter accordingly converges 
to zero. In order to reach this conclusion, thus, these theories require the space 
dimension to be less than or equal to 2. 

On the other hand, the method proposed by McBryan and Spencer [llj to prove 
decay properties of correlations in classical spin systems has been extended to bound 
the correlation functions in the one- and two-dimensional Hubbard models in f8], 
pro]. These theories make use of C/(l)-symmetry property of the system and deduce 
that the electron pairing-pairing correlation function for 2 separate sites x, y is 
bounded from above by 

„-Ci(e,-0y)+C2 ^ |tu.„|(cosh(e„-e„)-i) 

where Ci,C2 > are constants, {6'u} are arbitrary taken real parameters indexed 
by every site on the lattice, (tu,v)u,v is the hopping matrix and the sum with respect 
to u, V is taken over all the sites. A suitable choice of the parameters {^u} yields a 
decaying upper bound on the correlation function. As in [llj Koma and Tasaki [8] 
took {6'u} to be the fundamental solution of a Laplace equation on the lattice and 
concluded that the pairing-pairing correlation function decays as |x — y| — ;> +oo. 
Macris and Ruiz [T^ referred to a list of the possible parameters {^u} summarized 
in [T3] and extended the decay bounds obtained by Koma and Tasaki to be valid 
for the Hubbard models with long range hopping matrix as well. Since appropriate 
choices of {^u} have been found in one and two dimensions, this approach verifies 
the decay of the correlation functions in these low dimensional cases at present, to 
the best of the author's knowledge. 

Apart from these analysis to bound the correlation functions in low dimensions, 
Kubo and Kishi [9] proved that the susceptibilities for the Hubbard models at non- 
zero temperature are bounded from above by the inverse of modulus of the coupling 
constant in any space dimension under a few assumptions on the sign of the coupling 
constant and the lattice. 

In this paper we consider the equal-time correlation functions for a class of the 
Hubbard models at non-zero temperature and show that the correlation functions 
decay exponentially in the distance between the center of positions of the electrons 
and the center of positions of the holes in any space dimension if a norm of the inter- 
action term of the Hamiltonian is sufficiently small (see Theorem 12.41 and Theorem 
12.51 in Section [5] for the precise statements) . As in [5] , [TU] our approach essen- 
tially uses the [/(l)-symmetry property of the model. We start from characterizing 
the correlation functions as a limit of Grassmann integrals over finite Grassmann 
variables. The Grassmann integral formulations called the Schwinger functions are 
mathematical objects defined on a rigorous base. The [/(l)-symmetry property of 
the model is simply implemented in the Grassmann integral and enables us to con- 
vert the formulation into a form of multi-contour integral of the Schwinger function 
with respect to newly introduced complex variables contained in the covariance ma- 
trix. Our objective is, thus, set to find an upper bound on the Schwinger function 
inside the multi-contour integral. The evaluation of the Schwinger function is done 
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perturbatively. We find a volume-independent upper bound on each term of tlie 
Taylor expansion of the Schwinger function with respect to the interaction around 
zero. We require the interaction to be small so that the perturbation series of the 
Schwinger function converges. 

One advantage of this approach is that the space dimension of the system causes 
no technical difficulty and the resulting decay bounds on the correlation functions 
are valid in any dimension. However, as we need to go through the perturbation 
theory with respect to the interaction, the additional assumption is imposed on the 
magnitude of the interaction. 

Rigorous frameworks have been developed to control the perturbation theory 
for many-Fermion on lattice independently of the volume factor. In 1,5: Pedra and 
Salmhofer extended the notion of Gram's inequality to be applicable to the de- 
terminant of the covariance matrices appearing in many-Fermion systems. Their 
abstract theorem I15( Theorem 1.3] is general enough to cover the modified covari- 
ance coming into play in our construction and bounds its determinant independently 
of volume and temperature. The tree formula for partial derivatives of logarithm of 
the Grassmann Gaussian integral summarized in |17) coupled with the determinant 
bound of the covariance makes it feasible to find volume-independent upper bounds 
on the perturbative expansion of the Schwinger functions. Though in this paper 
we employ the tree formula to bound the Schwinger functions, the same goal can 
also be achieved via a concise representation of the Schwinger functions established 
by Feldman, Knorrer and Trubowitz in [5]. Indeed, their bound [21 Theorem 1.9] 
proved in a quite general context needs only the determinant bound and the L^- 
bound of the covariance, which are the same inputs to return the upper bounds on 
the Schwinger functions by means of the tree formula. 

We often refer to the recent article [7] for basic lemmas needed in our construc- 
tion. In fact, this work should be regarded as a continuation of which intended 
to explain mathematical tools to analyze the perturbation theory for many electrons 
in detail. 

As in [7] we directly treat the perturbation theory of the original model without 
introducing any multi-scale technique. Though the results can be presented explic- 
itly in a simple manner through the single scale analysis, it costs the temperature 
dependency of the interaction. In our analysis the norm of the interaction is re- 
stricted to be less than a constant times jS^'^^^ for the inverse temperature /3 in 
d-dimensional case. As a renormalization group analysis in the theoretical front let 
us remark that for a wide class of the Fermionic lattice models various mathematical 
properties of the correlation functions have been intensively studied by Pedra [13] 
in a generalized form. Pedra's multi-scale analysis remarkably concludes that the 
correlation functions can be qualitatively analyzed if a norm of the interaction is 
bounded by a constant times (log/3)~^ in a 2-dimensional case, or less than a con- 
stant times /3" with some n £ {—d — 1, —d/2 + 1/2] in d-dimensional cases (d > 2) 
as well. 

This paper is organized as follows. In Section [5] we define the model Hamilto- 
nians, prepare notations and state the main results of this paper. In Section [3] the 
correlation functions are formulated into a limit of the finite dimensional Grasss- 
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mann integrals. In Section 0] upper bounds on each term of the Taylor series of 
the Grassmann integral formulations are obtained. In Section [5] we first prove that 
the covariance matrix fulfills the necessary requirements for the Grassmann inte- 
gral formulation to be bounded perturbatively. By using the perturbative bounds 
on the Grassmann integral formulations, we then complete the proof of our main 
results on the exponential decay property of the correlation functions. Appendix Rl 
shows that the coefficient function defining the interaction term can be replaced by 
a unique anti-symmetric function. Appendix [B] presents a proof of existence of the 
thermodynamic limit of the correlation functions. 

2 Model Hamiltonians and main results 

In this section we define Hamiltonian operators and correlation functions together 
with notations used in this paper and state the main results. We use the standard 
terminology concerning the Fermionic Fock space without providing the definitions. 
They arc documented, e.g, in the book ^ or briefly in ^ Appendix A]. 

2.1 The Hamiltonian operator 

We are going to define our Hamiltonian operator on the Fermionic Fock space on the 
d-dimensional hyper-cubic lattice F — (Z/(LZ))'* (i, d g N) and the spin coordinate 
{t, 4.}. The lattice F* of momentum is given by F* := ((27rZ/L)/(27rZ))'^. We admit 
conventions that Kronecker's delta i5x.y takes 1 if the element x is identical with y 
in the set they belong to, otherwise and the function Ip of a proposition P returns 
1 if P is true, otherwise. For any vectors a = (ai, • ■ • , a„), b = • • • , 6„) of 
algebra, let (a, b) denote the sum ^et || ■ ||rii be the Euclidean norm in 

K", (•,•)(-„ be the inner product of C" and || • ||c" denote the norm of C" induced 
by (•, •)(-.„ . For any finite set S let US' denote the number of elements of S. 

With the creation operator -0*^ and the annihilation operator Tp^^ for any x G F 
and ^ G {t, i} the free part is defined by 

x,yere,0e{t,i} 

where the short range hopping matrix {r(x^, y0)}(x,5),(y,0)Grx{t,4.} is given by 
/ d 

- t' ■ ld>2 ^ ((^x.y-ej-efc + f^x.y-e^+efc + ^x^y+ej-e^ + (^x.y+ej+efc) ^ ^J■3x,y 1 • 
j.k=l ) 

j<k 

The real parameters t, t' and fi are called the nearest neighbor hopping amplitude, 
the next to nearest neighbor hopping amplitude, and the chemical potential, re- 
spectively. The vectors ej G F (j' e {1, • • ■ ,d}) are defined by ej{l) := Sjj for all 
jjG {!,••• 
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Since our problem becomes trivial otherwise (see Remark l2.10p . we assume that 

\t\ + \t'\U>2^0 (2.1) 

throughout the paper. 

To define the interacting part of the Hamiltonian we introduce functions Ui : 
{Z'y X {t,i}^' ^ C (/ G {1, • • • ,"}) satisfying the equality 



C/;((xi,--- ,x/),(,^i,-- - , (-2 2) 

= c/,((xi,... ,xo,(0i,-- - - ,6)) (v;>i) 

and the translation invariance 

C/;((xi,X2, • • • ,xO, (6,6, • • • (01>2, • • • ,0/)) 

= C//((xi+y,X2+y,--- ,x,+y),(a,6,--- , 6), ((^i, '^2, • • • >0) (VZ > 2) 

(2.3) 

for aU (xi,... ,xOe(Z'^)', (6, • • • , 6), • • • , 0/) e {t, 1}' and y e Z-^. 

We define a restriction of Ui by periodicity in the following way. Let [L/2J 
denote the largest integer not exceeding L/2. For any x e Z'' there uniquely exists 
x^ e ({-[L/2J,-[L/2J + I,-- - ,-L-/-/2j + L - l}y such that x x^ in T. By 
using this identification we define Ul.i '■ (Z"^')' x {^, J,}^' C (Z G {1, • • ■ , n.}) by 

;7l,i(x,6 0) :=C/i(x^,60), 

C/l,/((xi,X2, • • • ,X/_i,X/), (6, • • • ,C;)> •• • ,0/)) 

:=C/,(((xi -xO^,(x2-xO^,-- - ,(x,_i - xO^,0),(Ci,-- - ,CO,(</'i,-- - ,M) 

(2.4) 

for ? G {2, • • • ,n}. Note that [7^^/ is periodic with respect to the spacial variables 
and obeys ((2?2)) - ((2?3)) and 



^lim Ul.i{{xi,--- ,x,),(6,-- - ,6),(0i,-- - ,0/)) 
L^iT (2.5) 

= C/;((xi,--- ,x,),(6,-- - ,6),('/'ir-- ,00) 

for aU (xi,... ,xj) GZ', (^i, . . . , ^;), (,/,i, . . . , 0,) g (W G {1, ■ • ■ , n}). 
The interacting part V is defined as follows. 

n 

^^=E E E ,xO,(Ci,--- ,6),('/'i,--- --/-O) 

/=1 xj-er jj,.*3e{T,4} (2.6) 
vje{i, -- ,1} vje{i, .. 

Note that the condition (|2.2p makes V self-adjoint. The following examples motivate 
us to generalize the interacting part of the Hamiltonian as defined above. 
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Example 2.1 (the density-density interaction). With real functions Uf'^ : (Z'')' x 
{t; iY — > R (/ G {1, ■ • • , n}) satisfying the translation invariance 

t/,'^'^((xi,... ,xo,(a,-- - ,6)) = ;^i''((xi+y,-- - ,x, + y),(6,-- - 

and U^^^i^i,- ■ ■ ,xz), (^1, • • • ,6)) = if (x„e,) - (xfc,^) with j ^ k {yi > 2), let 
US define the density-density interaction Vdd by 

n I 

E E E • ■ • '^')' (^1' • • • '^0) n ^;.C.^x,c„ 

V3e{i, -- ,!} V3e{i, -. ,i} 

where the function Ul'^i is derived from Uf'^ in the same way as in (|2.4I) . The 
operator V^d can be rewritten as 

n i 

Vdd=Y. E c^l'((xi,--- ,xo,(ei,--- ,6))n'^c.A 

/=1 XjET ij,4'j€{t,l} j = l 

vje{i, -- ,1} vje{i, .. ,o 

which shows that Vdd has the form (|2.6p . 

Example 2.2 (the spin operator coupled with a local magnetic field). Introduce a 
local magnetic field = {B^J^\ B^^^ ) ■.Z'^^R^ and let Bl,^ = 
(^l L^ ^i^L' -^if !c) : r ^> be the restriction of Bx on T by periodicity as defined 
in (Hai). With the spin operator = {sP,S^\s^^) (x G T) given by 

S^^:=l E OxC'/^x^ (^e{l'2,3}) 



2 

«,0e{t,i} 



with the Pauli matrices 



we define a self-adjoint operator Vs by 

K :=^(Bi,x,Sx). 

xer 

Since 



and P^^'^ = P^'^ (VZ £ {I, 2, 3}, V^, (/) e {t, i}), the operator Vs provides one example 
of the interactions of the form 
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Example 2.3 (the spin-spin interaction). With the spin operator Sx introduced in 
Example l2.2l and a real function w(x) : Z'' — ^ K, we define the spin-spin interaction 
Vss by 

Vss ■■= ^ WL(x-y)(Sx,Sy), 

x,yGr 

where the coefficient wl is the restriction of w on F by periodicity. A calculation 
shows that 

Xl,X2Gr5l,52,01:02G{t,i} V ' = 1 / 



Hence l^^s can be written in the form (j2.6p . Consequently, our V covers the inter- 
action of the form V^d + V^s -l- Vss ■ 

In this paper we treat the Hamiltonian operator H = Hq + V. 
2.2 Main results 

We employ norms || • and || • ||/ to measure the magnitude of the interaction. 
\\Ul.i\\l,i ■=max max |C/l.i(x, ^, 0)|, 

max max E E E 

vfce{i,- - ,1-1} vkeii.- - vke{i.- - ,1} 

■|[/L,i((xi,... ,x,_i,0),(a,--- (V?>2), 

lC/i||i := sup niax V |C/i(x, ^, 

0e{t4} 



"^'ll'^=,,??^^ne.^T;} E E E 

•|C/z((xi,... ,x,_i,0),(a,-- - ,6),('/'i,-- - (W>2). 
For any operator O we define the thermal average (O)^ by 

Tr(e-'3^C') 



where the trace operation is taken over the Fermionic Fock space on F x {'I", J,} and 
the positive constant /3 is proportional to the inverse of temperature. 
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Define the function Ft^t',d '■ ^ ^>o by 



2{\t\ + 2{d-l)\t'\) y A{\t\ + 2{d-l)\t 



n\2 



+ 1. 



(2.7) 



For X G Z'' let us define tp^^ by considering x as a site of F by periodicity. 



The main results of this paper are stated as follows. 
Theorem 2.4. Assume that there exists R G (0, 1) such that 



l/(2e7i-d) 



+ 1 



l/(2e7rd) 



R. 



(2.8) 
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For any to S N and any Xj, G Z'', ^ji "Aj G {t) i} S {1, • • ■ , m}) f/ie thermo- 
dynamic limit 



-LSN 

exists and satisfies the following inequality. 



(2.9) 



lim 

-L-H-oc 



< (4™+i - m42™+i log(l - R)) ■ Ft,t'4 ( ^ 



(2.10) 



In the case that the interaction V is the on-site interaction 
^ Sxer V'xtV'xj.'^/'xJ-V'xt {U G M), the 4 point correlation functions can be bounded 
more sharply as follows. 

Theorem 2.5. Assume that the Hamiltonian operator H is given by 

H = Ho + UY, V'xtV'x^V'x^V'xt 
xer 

with the coupling constant U €M. satisfying 



\U\ < (108^)-^ 



l/(2e7rd) 



(2.11) 
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For any xi, X2, yi, y2 G the thermodynamic limit 

^lim (V'iit^isiV'ysiV'yit + ^U\^*y2i^^2i^^i^) L (2-12) 



+ + 



exists and satisfies 

-4i3ll'<^i+*2-yi-y2|lKd 



< 324 • Fi 



t,t',d 



Theorem [23] and Theorem 12 . 5 1 will be proved in Section [5j 
Remark 2.6. Theorem 12.51 does not follow Theorem 12.41 To see this, we write 



xer 



/c((xi,Ci),(x2,6),(yi,0i),(y2,</>2))V'*i^iV'x2e2^yi</'i^y202 



Vje{i,2} vje{i,2} 



with 



/c((xi,^l), (X2,6), (yi,0i), (y2,(/>2)) 

■= ^ J^xi,X2^yi,y2'5xi,yi('55i,t%,i ~ '5?i,i%:T)('^0i,t'^024 " S^i,iS^2,t)- 

The function /c satisfies the anti-symmetricity (jA.ll) and 
max-^ max V |/c((xi, ^1), (x2, 6), (yi, 0i), (y2, 02))1, 

vje{i,2} 

max V |/c((xi,^i), (x2,6), (yi,0i), (y2,02))| > 

(yi,0i erx{t,i} -"^ — ' 

vje{i,2} 

= M 

2 

Thus, Lemma [A. II proved in Appendix El ensures that if a function J7c : (Z^)^ x 
{t,!}'' ^ C satisfies (E^l-dlSl) and its restriction Uc,l ■ x {t,i}^ ^ C by 
periodicity obeys 

U Y V'it'^x^V'xi'f/'xt 

xer 

^ Y C^c,L((xi,X2),(Cl,6),('/'l>2))V'ii{iV';252V'x202V'xi0i, 

xi,x2erei,52.<?ii.<?i2e{t4} 

(2.13) 
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the inequalities 

^<I|c/c,lIIl,2<I1c/c||2 (2.14) 

must hold. 

Let us apply Theorem 12.41 to the interaction (|2.13p . The inequalities (j2.8p and 
(|2Ta imply 



\U\ < (256/3) 



l/{2eTTd) 

- 1 



which is a stricter constraint than (j2.11l) . Moreover, for R close to 1, the coefficient 
43 _ 2. 45 i^g^i _ larger than 324. 

Remark 2.7. One can prove that for any 6 > there exist constants Ct^t' ,d,b > 
depending only on t, t',d,b and ^, ^ > depending only on t, t' , d such that 



Ft.t' ,d ( 5^ I +1 



l/(2e7r<i) - ^"t,t',dP 



- 1 



for any P >b. Hence, the interaction needs to be small to claim the decay bounds 
in low temperatures. 

Remark 2.8. To generalize the results to many-Fermion systems with finite coor- 
dinates of colors is straightforward. We present the results only for the spins {t, i} 
in order to refer to proved materials for many-electron systems summarized in [7] . 

Remark 2.9. We use the translation invariance (|2.3p to prove the existence of 
the thermodynamic limit (j2.9D in Lemma 15.181 Without assuming (|2.3p we can 
also prove the inequality of the form (j2.10D with lim sup^_j._|_oo ^gpj in the left side 
instead of limL_j.+oo,LGN under an appropriate modification of the norm of Ui. 

Remark 2.10. If \t\ + |t'|lci>2 = 0, the correlation functions decay trivially. To 
prove this, let us take any real parameters {^xjxer and define the unitary operator 
Ab by 

Ag:= W e'''-'^;«'^='«. 
(x,c)erx{t,i} 
In this case, AgHAg — H and thus 

(2.15) 
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If LjLi % 7^ LfLi in r, we can ch oose {^xjxer to satisfy X^^i ^'x^ -LjLi %j ^ 
in R/(27rZ) and the equahty (12.15^ imphes that 

3 Grassmann integral formulation of the correla- 
tion functions 

In this section we formulate the correlation functions as a limit of finite dimensional 
Grassmann integrals. To attain this goal, we follow steps. As a preliminary let us 
fix the way to abbreviate the notations. 



3.1 The correlation functions 

To simplify presentations, from now we write = (xi, • • • , x;), = (yi, • ■ • , y;) G 
{Z'^y, = (^i,-- - ,6),*' = (01,- •• ,4>i) e {t,;}' (VZ G N). To indicate the sites 
and the spins on which the correlation functions are defined, we use the notation " 
and write 

S" = (li, • • • , U), l>™ - (01, • • • , 0™) e {t, ;}". 

We identify G (Z'*)' as an element of F' by periodicity without remarking 
when we are considering a problem defined on F'. 

To derive the correlation functions systematically, we introduce real parameters 
{A(X",y™,S™,$™)}_^„-, ;^„-,gp,i, 2^,^^g{^_^j„-, and define the cocflicient 

Ux,i{X\Y\^\ $') ■.^li<nUL,i{X\^','^')Sx,.Y, 

+ 1/=a(a(x*, r", s", + A(r", X", s")) 

for aU X',y' G F', G {t,4.}', ^ G {1, • • • ,max{m,n}}. We see that 

c/A,/(x',r',s',$o = c/A,i(y',x',$',s'). (3.1) 

Let us modify the interaction V to contain the coefficients {Ux^X^ ^ S', $')} 
and define 

max{m,n} 

'=1 (x',r',H',4>')er2'x{t.i}^' 
We set Hx Hq + V\, which is self-adjoint by the equality p. II) . Also note that 



Hx 



TT 
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Lemma 3.1. 



loe 



9A(X™, y™, S™, $™) \Tre"^^o 



(3.2) 



Proof. The proof is parallel to that of [Z, Lemnia 2.1], based on [71 Lemma 2.3]. □ 
3.2 The perturbation series 

In order to characterize the correlation functions as a limit of finite dimensional 
Grassmann integrals, we proceed in the following steps. Firstly we expand the par- 
tition function Tt e~^-^^ / Tr e~^^° into a perturbation series of the variables 
{[/a,;(^', Y\E}, <&')}. Secondly we replace the integrals over [0, /3) contained in the 
perturbation series by the Riemann sums to derive a fully discrete analog of the 
perturbation series. We then show that the discretized perturbation series converges 
to the original one by passing the parameter defining the Riemann sums to infinity. 
The discretized perturbation series is formulated into the Grassmann Gaussian inte- 
gral involving only finite Grassmann variables. Combining the Grassmann integral 
formulation of the discretized partition function with the equality p.2p completes 
the characterization. 

The first step results in the following proposition. 

Proposition 3.2. 

Tre-'^^^ _ ^ (-ir 



7n—l 



m I max{m,n} ^ 

k=iy (xl'',Yl>',s[>',<;>[>')er^'>'x{t,ir-''' ° 

det(C((x^s)p, (y(/>s)g))i<p,g<52I!Liifc' 



(3.3) 



where := (xfc,i, Xfc,2, • • • ,Xfc,iJ, Y^" := (yfc,i, yfc,2, • • • ,yfc,ij; 
'■= (Cfc,i,Cfc,2, • • • ,?fc,ij, := {<l>k,i,<l>k,2,- ■ ■ ,<Pk,iJ, and 



(x^s)p := {xy+i^u,^v+i,u:Sy+i), {y(j)s)p := {yy+i^u: 4'v+i,u, Sy+i) (3.4) 



for p = J2k=i^k + u, u € {!,■■■ ,lv+i}, u e {0, • • • ,m - 1}. Th 
C(x^x, yipy) is given by 



e covariance 



C(xex, y^y) 'jf (^i^ _ (3.5) 
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(V(x, ^jx), (y, (/),?/) G r X {tj-H X [0,/3)) with the dispersion relation 



d 



£;k := -2i^cos((k,ej-)) -4i' • ld>2 ^ cos((k, e^)) cos((k, efe)) - ^. (3.6) 

7 = 1 j,'!=l 

Proof. For any operator O defined on the Fermionic Fock space, let 0{s) denote 
gsHoQ^-sHa for s e M and (O)^^ denote Tr(e-/'^«0)/Tre-'3^o. We can apply [3 
Lemma B.3] to derive the equality 

^-PH, ^ ^-pH, ^ ^-m. y ^ ^_iyn I . . . ds„l,,>...>,„FA(si) • • • VxiSm), 



which leads to 

Tre-^^^ ^ /■ 

m— 1 



n 



max{m,n} 



• (V'xi,!?!.! • • • ^xi,,,Ci,ii (Sl)'/'yi,,i0i,,i (Sl) • • • V-yi, 101.1 (Sl) 

(3.7) 

By recalling the definition [71 Definition B.2] of the temperature-ordering operator 
and its properties [3 Lemma B.7, Lemma B.9], we can deduce p.3p from p.7p . The 
covariance C(x^x,y(/)y) and the dispersion relation E'k have been proved to have 
the forms dS^D-dSisi) in 7!. Lemma B.IO]. □ 

As the second step toward the Grassmann integral formulation, we introduce 
the discrete analog of the expansion p.3p . By taking a parameter h £ 2N//3, we 
define the discrete sets [0,f3)h and [—/?,/?);, by 

[0, /3). - {O, i J, • • • , /3 4} , [-/?, /3). -/3 + i, . . . , /5 - 1 

Note that tt[0,/3)/i = /3ft., — 2/3ft,. If the temperature variables x,y are 

confined in [0,j3)h, we write 

C?i(xCa;,y(/)y) = C(x^a;, y^y). 
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We then define {Tr e''^"'- / Ti e~'^"°)h by 
•=!+ > ^ T- 



m ( max{m,n} 

fc=i y ifc=i (x^'=,y^''=,H5^'=,*|^'')er2'fc xit^i}^'*" Sfce[o,^*)fiy 
■ det(Cft((x^s)p,(y0s)g))i<p,,<^^>^^,^, 

(3.8) 

where the variables (x^s)p, (y0s)p G F x {t, i} x [0, /3)/i are defined by the same rule 
as (j3.4[) . Since any determinant made up of the elements Cft(x^a;, y(/)y) vanishes if 
the size of the matrix exceeds 2L'^f3h{— jJF x {f, 4,} x [0, (3)h), the sum with respect 
to m is taken only up to m = 21"^ (ih m 



Remark 3.3. The diagonalization of the covariance matrix 

(Ch(x^a;,y(?!)y))(x,e,:E),(y^0,y)Grx{t,i}x[o,0)h was presented in [3 Appendix C] for any 
/i e 2N//3. To refer to this resuh we take h from 2N//3. 

The following lemma states that the partition function Tr e^*^^^/ Tr e~^^° in 
Lemma [5m can be replaced by (Tr e"'^^^/ Tr e^'^^o)/!. 

Lemma 3.4. For any r > there exists TVq e N smc/i i/iaf Re(Tr e^f^"^ / Tr e^l^"°)h 
> /or an e 2N//3 wif/i /i > 2No/l3 and all X{X"\ F™, S", $*) e M mt/i 
|A(X™,r™,S™,$™)| < r. Moreo7;er, 



(?/'! c • • • V'! ; -0- I • • • 'Z'- I + ■0- I • • • "0- i "0- ; • • • 0- ; ) 



L 



1. a , /Tre-'3^^ 

= — — lim — — — — — — log 



B fe^ + =o f)\(Xrh yrh ■^m V Tr e"'^^'' 



(3.9) 



where for z G C wii/i Re z > 0, log z is defined by taking the principal value; 
log z := log |z| + z Arg z, Arg z G (— 7r/2, 7r/2). 

Proo/. By considering {C/a,;(^', S', $0}iG{i, . ,max{m,ri}},(x',Y',H',<E.')6r2' x{t,i}^' 
as mutually independent, complex multi- variables, we define the functions 
P{{Ux,i{X\Y\E\¥)}) and Phi{Ux,i{X\Y\E\<S?^)}) by the right hand side of 
p.3p and p.Sp . respectively. We show that P/i converges to P as ft. ^> +oo locally 
uniformly with respect to the variables. 

Pedra-Salmhofer's determinant bound [HI Theorem 2.4] (see also Proposition 
15.31 below for an extended statement) implies that 

I det{C{xp^pSp,yq(j)qtg)) < 4" (3.10) 
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for any (xp, ^p, Sp), (y^, 0p, ip) e T x {t,|} x [0,^) (Vp G {1, • • ' ,n}). Thus, if all 
the variables satisfy the inequality \U\,i{X'', Y\ S', $')| < r for an r > 0, we have 



(-1)' 



m / max{m,n} 

n E 



E 



I /-p 

• n / difcdet(C((x^t)p,(y0t),))i<p,,<^-^^i^ 

"11 77 E det(Cft((x^s)p,(y</>s)g))i<p.g<2-™^^i 



max{m.n} 



fe=i 



(3.11) 

where (x^t)p, (yt/)^),, (x^s)p, (y^s), are defined in the same way as in (13.41) . The 
right hand side of (13. lip is sunimable over to G N. 

Fix any ) G r^'*" X {t,;}^'" (Vfc G {l,--- ,™})- We define a 

function g,, : [0, Z?)™ -J- C by 

where Sj G [0, 13)h satisfies tj G [sj, + 1/h) (Vj G {1, • • • , m}) and (x^s)p, (y(/>s)g 
(p, g G {1, • • • , X^feLi ^fc}) S'i's defined by the same rule as (13.41) . 

Since C{yi^s,y(j)t) is continuous with respect to the variables {s,t) G [0,/3)'^ 
almost everywhere, so is det(C((xft)p, (y0i)g))i<p,g<x]"L '^ith respect to 
(ti,--- ,i™)G [0,/3)™, and thus 

^^m^ ghih, ■■■ , tm) = det(C((x^i)p, (y0i)«))i<p,g<Er=i 'fc 

for a.e. (ii,-- - ,<m) G [0,/?)™. Hence, by p.lOp and the dominated convergence 
theorem for L^([0,/3)™), we have 



lim 



— lim 
= 0. 



m 1-13 

n / difcdet(C((xet)p,(y#),))i<p^,<5^.._^,, 

m ^ 

"11);; E det(Ch((x^s)p,(y<?!)s)q))i<p,5<52;!Li 

fc=l SfcG[0,,3)fe 

, 



(3.12) 
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By p. lip and p.l2p we can again apply the dominated convergence theorem 
for ;i(N) to show that 



Km sup 



I y^l '^l if)l 



)\<r 



Phi{UxAX\Y\E',<S>')}) 



(3.13) 



for ah r > 0. 

Since the multi- variable function P — is entirely analytic, the convergence 
property (|3.13p and Cauchy's integral formula prove that 

lim sup 

with \Ux ; (X' , ,<i>^)| <!- 



d 



P {{UxAX\y\~\<^')}) (3.14) 



9^7A,fe(x^r^s^$'■o ^ 

~PH{{Ux.i{x'y,^'M)]) 

for aU r > and any /fc G {1, • • • ,max{m,n}}, (X'=, F'^, S*^, $'=) G T^^ x {t,!}^''. 

Since Tre~''^^/Tre~'^^'' > 0, the uniform convergence property (I3.13P verifies 
the first statement of the lemma on Re(Tre~^^^/Tre~'^^f);i. 

Note the equality that 



5 



log 



aA(X™,f™,S™,$™) VTre-^^oy^ 



Ph{{Ux4X\Y',E',<i>')}) 



PHi{UxAX\Y',E',<i>')}) 



(3.15) 



Combining the convergence properties p.l3p - p.l4p and the equality p. 151) with 
(I32D yields the equality dSl]). □ 



3.3 The Grassmann integral formulation 

As the third step we formulate (Tr e~^-^^ / Ti e~^^°)h into a Grassmann Gaussian 
integral on the finite dimensional Grassmann algebra 

{''/'xfa;! V'x5a:}(x.5,2;)erx{t,i}x[o,^i)h- Basic properties of the finite dimensional Grass- 
mann integral have been summarized in the books [1], |16j . We assume that each 
element (x, ^,x) G F x {|, 4} x [0,f3)h is numbered so that we can write 

r X {t, ;} X [0, p)h = {(x^a;), | j G {1, • • • , N}} 
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with N := 2L'^l3h. _ _ 

The Grassmann integral J •(iV'(xf2;)^. ■ • ■ '^V-'(x?2;)i'^V-'(x^rr)iv ' ' ' '^'^{yLix)i is a hnear 
functional on the complex linear space C[V'(x{a;) ' i^{xix)j I J ^ {1, • • ■ > ^}] of Grass- 
mann monomials and satisfies that 



• #(x|a:)N • ■ • #(x«x)i#(x|x)n • • ■ #(x$x)i = 1> 
y ■^(x«x)j, • • •■^^(x«x)j, V'CxIxjpi • • • V'(x«x)p, 

• #(x|x)n • ■ • #(x«x)i#(x|x)n • • ■ #(x«x)i = 0- 



(3.16) 



if fc < TV or g < TV. _ 

For simplicity let V'x) V'x denote the vectors of the Grassmann variables 
(V'(xea;)ir • • ,'0(xCx)„): (V'(xCx)i , " ' " ,V'(xCx)a:), respectively. We also write dtpx = 
#(xCx)« • • • ^^^ix)i and dtpx = rfV'(x4x)„ • ' • #(xCx)i- 

For any /(V'x; V'x) € C[V'(x{x), > V'(x|x),- I j e {1, • • ■ , iV}] the value of the Grass- 
mann integral / f{tpx,'>px)dipxd4>x can be computed by linearity and the anti- 
commutation relations of the Grassmann variables. 

For any N x N matrix Gh with det Gh ^ 0, let Gh denote the 2N x 2N skew 
symmetric matrix 

f Gh\ 
\-Gi )■ 

Prom the definition and the assumption that h G2N/I3 we see that 

= (-l)^(^-i)/2(detG,,)-i = (detG,)-^ 

Definition 3.5. For any N x N matrix Gh with det Gh'^O the Grassmann Gaus- 
sian integral / -d/iGh {'4'x,'>Px) : , V'(x«x)j | j G {1, • ■ • , AT}] C is defined 
by 



/ 



f{i'x,^x)di^Gdi'x,'>Px) ■= - 



Je-ii(^x,M',G^'{^x,M')d^xd^ 



X 



Transforming (Tr e ''^^/ Tre ^^°)h into a Grassmann Gaussian integral essen- 
tially relies on the following equality. 



Lemma 3.6. 



/ 



^(xjx),-^ ■ • ■^(x«x),i^(x|x)„ • ■■i^{^$x)^J^^GH{^x,i^x) 

= det{Gh{{x^x)j^,{x^x)pJ)i<u,v<k- 
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Proof. This equality follows [T, Problem 1.13] by replacing the notations ipx, tpx 
by "0X7 i'x, respectively. □ 



Remark 3.7. Our definition of the Grassmann Gaussian integral Definition 13.51 
differs from that summarized in (4j Problem 1.13] or its follower Jj Section 3.2] 
and corresponds to the statement derived by changing ipx for ipx ^^'^ '0x for V'x 
respectively in Problem 1.13] and T, Section 3.2]. In \T, Proposition 3.7] the dis- 
cretized partition function of the Hubbard model was formulated in a Grassmann 
Gaussian integral. The formulation [7J Proposition 3.7] required the symmetry as- 
sumption on the coefhcients ?7x,y,z,w (x, y, z,w G P), namely C/x,y,z,w = C^z,w,x,y 
Under Definition 13. 5[ however, we do not need any additional assumption on the 
coefRcients U\j{X'',Y\'E},^^) to complete the desired formulation below. 

Lemma 3.8. 



h 



(3.17) 

where 

'■~ ^ ^xi5ia;V'x2C22; ' ■ ' ^x,4;a;V^y!0!a;^y!-i(/-i-ia: • • • V'yii/iia;- 

Proof. First note that det Ch ^ for any h e 2N//3 by ff! Proposition C.7]. By 
applying LemmaEiHto det (C/i((x^s)p, {y4's)q))i<p,q<j2'^^iik '^^ and using the 
fact that J Idfic^iipxy^x) = 1, we have 

^ -1+ E 



Xj.g-/3ffo 



h 1 

"' m.— 



m I max{m,n} 

n E E E 

fc=i y ifc=i (x^'=,y^''=,H5^'=,*|^'')er2'fc x{t,i}^''' '>ke[Q.f))hj 

• V'yi. 101,1^1 • ■ • V'yi.ii0i.,isi • • • V'y^.i^^.is™ • ' ' ^y,r.,i^<t>^,i^Sr^d^,cM X^^x) 
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\ m=l 
m / max{m,n} 

■H E E E 

fe=i V ifc=i (x^'',Y^'=,E:J^'',*5^'')er2'fcx{t,4.}^''' s^e[o,p)h 

which is ((XT71) . □ 

Here we complete our Grassmann integral formulation of the correlation func- 
tion. 

Proposition 3.9. 

\^xiji ^Xriij,ft ^ ymi/'rf, ^yi0i ^yi0i ^ym^m^X^i^ri, ^xiCi/L 



/ie2N//3 



e 



Ei=i E(. 



(3.18) 

Proof. The equality p.lSp is obtained by substituting p.l7p into the right hand 
side of p.9p and differentiating the Grassmann polynomial 

by the variable A(X™, F™, S™, $™) inside the Grassmann integral (see 

[3J Problem 1.3] for differentiation of Grassmann polynomials). □ 

4 Perturbative bounds on the Grassmann integral 
formulation 

In this section we find upper bounds on each term of a perturbative expansion of the 
Grassmann integrals of the form (I3.18P for fixed h G 2N//?. Here we generalize the 
problem. The covariance is assumed to be any matrix G/i : (Fx {t, 4,}x [0, /3);i)^ — ?► C 
with detGh ^ 0. 
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Set an integral of Gh by 

I? := max < max — \Gh(x(,x,Yd>y)\, 

I (yAy)erx{t,i}x[o,/J). /i ,,,, 

max — \GhiY(t>V,x£.x)\ >. 

(y,0,«)erx{t,i}x[o,/J). /i J 

(x,?,x)Grx{t,i}x[0,^)h 

We define the Schwinger funetion by 

(4.1) 

for any TO e N and (X™, f", S™, 4") € T^™ x {t, if"' and any r; G C for which 
the denominator of (j4.1l) does not vanish. In fact, since 



Sj^A yrh arh ^rh {G h , Tj) IS analytic with respect to in a neighborhood of origin. 
The purpose of this section is to prove the following proposition. 

Proposition 4.1. Assume that there exists a positive constant B > such that for 
any p e N, (x^, , x^J, (y^, , 0fe, , y^J S T x {f,!} x [0,/3)/i (V^ £ {I,-- - ,p]) and 

sup |det((u,,v™)c„G/i(xj,^j,Xj,,yfc„(?!)fc„yfe„)) < I <BP. 

ll"!llc".IK!llc"<i vie{i,... ,p} 

(4.2) 

Let hm denote the coefficient of rj"^ in the Taylor series of the function rj i— > 
Sj^rh ~rh j,rh {G h , Tj) around G C. The following hounds hold. 

|6o| <B"\ 

%n\ < ^ f ^M'B'-l||[/L,/||L,iPj (Vto G N). 
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Proof. Let us consider {?7A,i(^% ^\ ^'0};e{ir - ,max{m,n}},(x',Y',H',<E.')er2ix{t,;r' 
as mutually independent, complex multi- variables. In a neighborhood of origin we 
can define a multi-variable analytic function W{{U\^i{X'- ,Y\E'- ,^^)}) by 

W{{Ux,i{X',Y\E\^')}) 
Its Taylor expansion becomes 



„ ITT'- . 
TO=0 fe=l 



m 



m / max{7?i,n} 

E 

^ (x^^Y■,'^4^^.^'=)e^^'-=x{t,-^}'''=/ 

V/G{1,--- ,max{m,n}} 

v(x' ,h' )er^' X {t,4.}^' 



We can take a small r > so that for all i) € C\{0} with I77I < r 

S j^rh ■f'm ^ihif^ h •} 

1 9 



Z^'? l9J7A,m(^", S™, 



m —I I h 



\ 



E (m _ IV II I E E 



(4.3) 



n 



W{{UxAX\Y',E',^')}) 



V/£{1,'" ,max{7h,n}} 



fe=i 
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By comparing with the definition (j4.ip we see that the equahty (j4.3p also holds true 
for ?7 = 0. 

We will bound the coefficient of 77™ in the right hand side of (|4.3p . To organize 
the argument we divide the rest of the proof by 4 steps. 
(Step 1) First remark that by Lemma 13.61 

bo = 



-Wi{Ux,iiX\Y',E',¥)}) 



t/A,i(X',y,H',*')=0 
V/G{1,--- ,max{m,n}} 
V(X',F',H',$')er^'x{t,^}^ 



I3h 

det(Gh(xjCjO,yfe0/cO))i<j,/c< 



(4.4) 



Thus, by the assumption g2|) |6o| < 

(Step 2) Consider the coefficient of 77™ in (|4.3p for m > 1. Here we ap- 
ply the tree formula |17[ Theorem 3] to characterize the partial derivatives of 
W{{Ux^iiX\Y'- evaluated at origin. We then obtain 



ni I n 



/3 m! n I 



E 



d 



d 



n 



V/G{1,--- ,max{m,n}} 



m I n 



5] f/i,,(4^si^$L'= 

'-)er 



, A(M(T,7r,s)) 



•^0 



7rGS„ + i(T) 



fe=l 



X^*" .X^*" ' ^^-^ 



^''x = 't''x=° 

Vge{0,l, -- ,m} 



(4.5) 



Definitions of the newly introduced notations are in order; T({0, 1, • • • ,™}) is the 
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set of all trees on the vertices {0, 1, • • ■ , m}, S„j+i(T) is a T-dependent set of per- 
mutations over m + 1 elements, x(T, tt, s) is a (T, tt, s)-dependent non-negative real 
function obeying 

f ds x(T,7r,s) = l, (4.6) 

M(T, TT, s) is a (T, tt, s)-dependent non-negative real symmetric (m + 1) x (m + 1) 
matrix satisfying M(T, 7r,s)p_p = 1 for all p € {0, , m} and the Laplacian 
operators Ap^g and A(M(T, 7r,s)) are defined by 



\ * . . * 



Ap,, :=-((-^j (^^) ), A(M(r,7r,s)):= ^ M(r, tt, s)p,,Ap. 



9' 

P,9=0 



where 

d 



are vectors of the Grassmann left derivatives corresponding to the labeled Grass- 
mann variables {■0(xCx), - '^U^x)^ \ J ^ {^r ■ ■ , N}} (Vq G {0, 1, • ■ • , m}). 

(Step 3) We will bound the right hand side of (j4.5p by replacing the sum over 
trees by the sum over incidence numbers in the next step. As the preparation let us 
summarize some facts. Let do, di, • • • , dm denote the incidence numbers of a tree T 
corresponding to the vertices 0, 1, • • • , m, respectively. Every term of the expansion 
of J^Ip ^y^j^{Ap^q+Aq_p) IS & product of m Laplacians as the tree T has m lines. Each 
product of the m Laplacians contains dq derivatives with respect to the Grassmann 
variables with label q for all q G {0, 1, • • • , m}. Moreover, the dq derivatives consist 
of aq derivatives with respect to the variables {'0x^2,}(x.5.2;)erx{t.4,}x[o,^)h ^-nd dq 

derivatives with respect to the variables {i^i^x}{x,^,x)£rx{^,i}xlo,i3)h ^'^^ some non- 
negative integers aq, dq satisfying dq = aq + dq. 
This product of the m Laplacians produces 







' Ik 










] dk 






V <3fc . 





monomials when the term acts on the monomial 

m 

• n (v^i.c.. 

Sfc ^Xfc_24fc,2Sfc ' ' ' ^Xk,l^^k,lf.Sk (4-7) 

fe = l 

' ^yfc,ifc0fc,ifeSfc''^yfc,ifc-i0fc,ifc-iSfc ' ' ' ''^yk,i4'k,iSk^ ■ 

Moreover, every remaining monomial left after the product of the m Laplacians 
acting on (|4.7p consists of m — ag -I- X]fcLi(^fc ~ ctfe) products of the Grassmann 
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variables of the form V'x^s and rh — uq + ~ Q^fe) products of the Grass- 

mann variables of the form i/'x^s a^^d is acted by the operator e^'^^^^"'"'^'**-'''. The 
determinant bound (j4.2p and the non-negative symmetric property of M(T, tt, s) 
together with the fact that all the diagonal elements of M{T, tt, s) are 1 validate 
the inequality 

A(M(T,7r,s))-T0 _-tO ,0 ,0 



fc=l 



(4.8) 



vqe{o,i.- - .™} 



Z^g = "Q~^q = 

(see [71 Lemma 4.5] for a detailed proof of the same bound in the case B = 4). 

The coefficient of each remaining monomial left after the product of the m 
Laplacians acting on the monomial (|4.7p is a product of m elements of the covariance 
(Gft(x^a:,y(/)y))(x,^,2,),(y,0,2/)6rx{t,4.}x[o,/3);, and going to be multiplied by 
UT=i \ Ul,i, {Xl" , , ^[^ ) and summed over the variables (Xl" , ,'S>[\sk)€ 
r'*" X {t, 1}^''° X [0, f3)h for all fc G {1, • • • , m}. By considering that the tree T starts 
from the vertex and using the properties (j2.2p - (l2.3p . a recursive argument through 
all the lines of T running from the later generation to the earlier generation shows 
that this sum is bounded by 

m 

n iit^i.'ju.'.^™- (4-9) 

k=l 

Note that to derive the upper bound (|4.9p we repeatedly use the inequality of the 
form 

^ E E |C/L,K^',S',<f')||G„(yay,x,C,s)| < ||C/i,niL,;2?, 

se[o,/3)h (x',H',*')er'x{t.i}^' 

where (y,CT, y) e F x {t, 4-} x [0,f3)h is any fixed element and Xj and £_j are j-th 
component of X'- and 5' respectively for some j € {!,■ ■ ■ , 

Lastly remark that the number of terms containing a, derivatives with respect 
to the variables {V'xjx}(x,?,x)erx{t,i}x[o,/3)ft (Vg £ {0, ,™}) in the expansion 
n{p.g}eT(^P,9 + ^9,p) is at most 



n 

q=0 



dq 



(Step 4) By summing up all these considerations, replacing the sum over trees 
by the sum over possible incidence numbers da,-- - ,d„i satisfying the condition 
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X]^o '^9 ^ ^"^i ^^^'^ using Cayley's theorem on the number of trees with fixed 
incidence numbers (see, e.g, [22l Corollary 2.2.4]) and the equality (14. 6p . we have 



1 6. 



<-- y n 

soG[0,/3)h fc=l 



( ' 



E 



sup 

se[o,i]" 



TeT({0,l:-,ni}) -es„_^i(T) 



,A(M(T,7r,s)) 



Jl (Ap,, + A,,p) • (the monomial K7\ ' 



{p,9}eT 



V<je{0,l, -- ,m} 



1 1 



m / n 



(to - f)! 



soe[0,/3)h fc=l Mfe=l 



dfcG{l,--- ,2lk} 
Vfce{l,--- ,m} 



V TTi w f^™ ^™ ~ TT ' U3™-"o+i:r=i('fc-"'=) 



Vqe{0,l, -- 



9=0 



/ TO 




' TO 












V "0 









fc=i / 



fe=i 

1 ifc-m 











) "^-i 1 









E ^T,'^=0 "9=^7=0 "1 = 



afc! 



c«Oi|5oG{0,--- ,m} 
VfeGfl.--- ,m} 



9=0 



m \ m 



n 

fc=i 



h \ f I 



ctk J \ ak 



Then, by dropping the constraint ^^"^_g aq=J2'"-o ^^'^ using the equality 

I I 



EE("+") 



a=OQ=0 



I \ f I 



a \ a 



= 14^, 
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we obtain 



7n / n 



^ k=l 



m 

k=l \lt = l / k=l 



\i=i / 

□ 

The evaluation of the tree expansion presented above overcounts the combinato- 
rial factor, which is defined as the number of monomials appearing in the expansion 
Yl^p g}eT(^p.9 + ^g,p) ' (the monomial (|4.7p ). If we consider the 4 point correlation 
functions for the on-site interaction V = U X^xgf ^xf^xi^'x^V'xt {U G R), the exact 
calculation of the combinatorial factor without overcounting is possible as proved 
in [7J Lemma 4.7]. Consequently, the perturbative bound is improved in this case. 

Proposition 4.2. Assume that a non-singular matrix Gu satisfies (|4.2[) . For any 

m G NU {0} let Cm denote the coefficient of rj™ in the Taylor series of the function 

ri 1 — > 

around G C. The following bound holds. For any m G N U {0} 

\c^\<^f''''^')iVB\U\r. 
6m + 4 \ m J 

Proof. The argument in [7i Section 4] leading to the bound [71 Lemma 4.8] straight- 
forwardly applies to deduce the claimed bound. □ 



5 Exponential decay of the correlation functions 

In this section we prove the exponential decay property of the correlation functions. 
The Grassmann integral formulation p. 181) serves the base of our analysis. In 
the subsection 15.31 we will see that the J7(l)-invariance property of the Grassmann 
integral leads to a replacement of the covariance Ch by a modified covariance in 
the Grassmann integral formulation. In the following subsections 15.1115.21 we study 
properties of such perturbed covariance matrices as preliminaries for proving our 
main theorems. 
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5.1 Determinant bound on modified covariance matrices 



Here we prove the determinant bound on a generalized covariance matrix of the 
form 

Cs{^^^,yM := ^ E e'<---)e-(-^)-^ (if^ - (5.1) 

(V(x,C,x),(y,(/>,2/) e r X {t,;} X [0,/3)). In ^J^i fk : T* ^ C is any function 
satisfying 

llmfkl < ^ (Vker*). (5.2) 

Note that the inequahty ([Ol) imphes |1 + e^^"], |1 + e"^^"] > (Vk S T*), and 
thus Cg is well-defined. The covariance (j5.ip is a generalization of the covariance 
p.5p in which the dispersion relation £^k is real-valued. 

The volume- and temperature-independent determinant bound on the covariance 
(|3.5p was established by Pedra and Salmhofer in [TS] by extending the notion of the 
Gram bound on determinant. However, the Gram representation of the covariance 
presented in [TS] Section 4.1] does not apply to the case that the dispersion relation 
is complex- valued. We need to construct a Gram representation of the covariance 
matrix (|5.ip to bound its determinant. 

We use the following abstract framework proved by Pedra and Salmhofer as a 
generalization of Gram's inequality. 

Lemma 5.1. JTH Theorem 1.3] Assume that K,k e NU {0} satisfy k + K > 1. Let 
X be any set. Assume that J is a set totally ordered under H ' and is a strict 
total order in J. Let Q, be maps from X to J (V/ G {1, • • • ,k + K}). Let % be a 
Hilbert space with the inner product and || • W-u denote the norm induced by 

{■, ■)^. Define a matrix M^^y {x, y EH) by 

k K+k 

;=i i=k+i 

where f^^, gl^eTL. If 

supmax{||4||„,||g^||„}<7, (VZ e {0, • • • , + fc}), 

a:eX 

the following bound holds. For any m, n £ N and xi, - • ■ , a;„, i/i, • • • , j/„ G X 



2n 

sup |det((uj,Vfc)c™A/4^,yJi<j-fc<„| < ( ^ 7i ) . (5.4) 

||Uj||c™,||Vj||c,Ti<l 

VjG{l, ■■■,«} 



/k+K \ 



Remark 5.2. The determinant bound [151 Theorem 1.3] was originally claimed for 
the case that m — n in (j5.4p . The bound for general m immediately follows from 
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(|5.4p without the coefficient (uj,Vfc)j.„. In fact it is sufficient to redefine H by 
C™ (g) H and M^^y by 

1=1 

K+k 

for any vectors u^^'Vx £ C™ with ||u2.||c™, Hv^^Hc" < 1- 

As an application of Lemma l5.ll we prove the following statement. 
Proposition 5.3. // the inequality 

|im£:k|<^ (Vker*) (5.5) 

holds, then for any m,n and {xj,^j,Xj), {yj, (pj, yj) £ F x {t, 4-} x [0, l3) (Vj e 
{!,••• 

sup I det((uj, Vfc)c™Cf (xj-^jXj,yfc(?!)fc?/fc))i<j,fc<„| < 4". (5.6) 

l|ujllc-,l|vj||c"<l 
ViG{l,--- ,«} 

Remark 5.4. The upper bound 4" claimed above is same as that obtained in [TSJ 
Theorem 2.4] for the covariance matrices with real- valued dispersion relations. 

Proof of Provosition \5.Sl Our argument here closely follows ^15t Section 4.1]. Since 
r* is a finite set, for a sufficiently small e > we may assume that Re(fk + £) 7^ 
for all k G r*. Set f k - ^^k + e (Vk e F*) and define the parameterized covariance 
CI by 

Sr ,i(k,y-x)^-(y-x)ek,e f h-x<0 _ h-x>0 



Cli^^xMy) ■■^ -TtT e'^<k,y-x) tlZl^ 



ker 



We construct a Gram representation of C| fitting in the form (j5.3p to apply Lemma 
5.11 Observe the following decomposition. 



C|(xCa;,y0j/) 

- k—>OJJ e (^lRo£k,.>0^^g_^£^^ + lRc£k..<0 ^ ^ g/3£^,, j 



(5.7) 
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Note that for i > and f G C with | lm£\ < vr//? and Rc£ > 



-e 



-tReS 



l + e-f^£ |l + e-'3^|2 

(5.8) 

|l + e-/5^I2 ''t-2 + (Re£:)2' 

where f denotes the complex conjugate of £. 

By substituting the formula (|5.8p into (|5.7p we obtain 

ker* 

• ly-x<0 1rc£,,.>0 h , ,-/3£.,.|2 ^ / '^^I 



/ e-'(2'-^)i'"^'-.-(l + e-^^) Refk,^ /" , e*^^^"^)" 

y |l + e P^-k,e|2 7r u2 + (Re£:k,e) 

g~i{l3+x-y)I-m(-eu,e) (^l _|_ g/3£k,£^ Re( — £^k e) f g«(/3+2;-i/)u 



\2 



(5.9) 

Let us define the Hilbert space L^{r* x {t,i} x K) equipped with the inner 
product (•,-)i2(r.x{t,i}xK) given by 

(/'.9)L2(r*x{t4}xR) - I dvf{k,^,v)g{\<i,^,v). 

For (x,^,x) G r X X and a G {-1,1} define vectors f^^^, ^ G 

L2(r* x{t,;}xM)by 



l^^^fSaSy,.. /I Refuel 



iRe(a£,,.)>oO5,0e 1 1 + g-ZSaf.,, | V ^ i^, + Re(a£k,s) ' 



|1 _,_ g~/3a£k,.| V w + Re(a£:k,£) ' 

Note that 

(/x,{,£!;'/y,0,y)L2(r*x{t.i}xR) = (/x,C,a;'5y,0,y)L2(r*x{t4}xR) .^^^^ 

= (5x,5,a:;3y,0,y)L2(r'x{t:i}xR) ~ 0, 
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and 

ll/x,C,a;llL2(r'x{t.i}xR) = ^ X/ -'-Ro(a£k.. )>0 1 

_ 1 1 ^'-''^ 

\\9i,iJ\h{r'x{t.l}xR) - JlYl lRo(a£k,,e)>0 _^ g-/3a£:k,,|2- 

The equalities (ICTU)) - (l5lT|) imply that 

11/^,4,2: + /y,0,j^lli2(r-x{t,i}xR) ^ ]^ X! (lRc(oek,e)>0 + lRc(-a£k,e)>o) = 1 (5-12) 

and under the assumption that | Imfkl < / i'^P) 
Il5x,c,x + 5y,0,j;lli2(r* x{t4}xR) 
< 



ker* \ 



1 

>0 



g-2/3Re(afk,e) + 2 C0s(/3 Im £:k,e)e"'^ + 1 

1 



+ lRc(-a£k,e)>Og2/3Rc(a£k,.) + 2 COs(^ Im f k,e)e'5 I^'=('^^>'-) + 1 ^ 

1 / 1 1 \ 

- Jd I lRo(a£k,e)>0g_2,3Rc(a£k,.) -|- 1 -'"R<=(-"£k,s)>0 ^2/3 Ro(a£k, J _|_ 1 ) " 

(5.13) 



ker* 



for all (x, e, x), (y, y) G T x {t, 1} x [-/3, /3] and a e {-1, 1}. 
The equality (15.91) can be written as 

C|(xfa;,y0y) = la-i:<o(/x,^, 1,-/3 + /xi,-x' 5y,</,,y + 5y,0,-y)L2(r-x{t,i}xR) 

+ ly-.T>o( ^ /x.?..T ^ /x,4,-K'5y,0,a + 3y.0./3-;^)L2(r*x{t,i}xR)- 

(5.14) 

The representation (j5.14l) and the bounds (|5.12p - (|5.13p enable us to apply Lemma 
15.11 to obtain the bound (|5.6p for C|. Then, by sending e \ and the continuity 
of C| with respect to e we complete the proof. □ 



5.2 Exponential decay of the covariance 

In this subsection we show the exponential decay property of the covariance and 
find a volume- independent upper bound on its integral as the corollary. 

Later in our analysis we will need to deal with a modified covariance whose 
dispersion relation is given by -Bk+ze with z £ C, p G {1, • • • ,d}. It is convenient to 
prove the L^-bound on the modified covariances which include the normal covariance 
(|3.5p with the real-valued dispersion relation p.6p as a special case at this stage. 
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Furthermore, to bound the covariance with the dispersion relation iJk+zep needs 
to control the covariance with the perturbed dispersion relation i?k+ze +we defined 
by 



Ek+zep+wBg ■= -2t ^ cos((k + zBp + weq, ej)) 
d 

-W ■ld>2 ^ cos((k + zGp + weg,ej)) cos((k + + wGq, Gfc)) - ^ 



i,k = l 
j<k 

for z,w € C and p, g G {1, • • • , d}. Let us study properties of -Ek+zep+uie, • 
First we clarify a sufficient condition for the inequalities 

7T TT 
Im(i?k+2ep+u)e,)| < ^, I Im(i?k+zep+tue,)| < 

which will ensure the analyticity of the modified covariance with respect to the 
variables z,w and make the determinant bound Proposition 15.31 applicable. 

Lemma 5.5. Let z,w G C and r > 0. // | lmz|, | Imw| < ^ log Ft ^t',d{r) with the 
function Ft^t',d{') defined in (12. 7p . then | Im_Ek+zep+iue, | < f holds for all k G F* 
and p^q G {1, • " ' id}. The parallel statement with '<'s in place of '<'s holds as 
well. 

Proof. Let us write z = a + bi, w = u + vi (a, b,u,v eM). By using the formula 

cos(a; + iy) = cos{x) cosh(?/) — i sin(x) sinh(?/) {x, y G M) 
and writing k ~ (fci, • • • , kd), we find that 

d 

ImSk+zep+iwe, = sin(fcj + aSpj + uSgj) smh{bSp^j + vSqj) 

i=i 

d 

+ 4t'ld>2 ^ cos{kj + aSpj + uSqj) cosh{b5pj + v6qj) 
■ sm{ki + aSp^i + u5q,i) sinh(Mp,; + vSq^i). 



j.i=i 



Assume that \v\ < s for s > 0. 

I Imi?k+zep+u)e, I 

d d 

< 2|t| ^sinh(s((5pj + Sqj)) + 4|t'|ld>2 ^ cosh{s{Spj + (5gj)) sinh(s((5p,/ +(5g,;)) 

1=1 j.i=i 

< lp=q{2\t\ sinh(2s) + 4|t'|(d - 1) sinh(2s)) 

+ lp5^g(4|i| sinh(s) + 8\t'\{d - 1) cosh(s) sinh(s)) 

< 2sinh(2s)(|t| +2(d- l)|i'|), 
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where we used the facts that for s > 

1 < cosh(s), 2 sinh(s) cosh(s) = sinh(2s), 2sinh(s) < sinh(2s). 
Thus, the inequahty 

2sinh(2s)(|i| +2(rf- l)|i'|) < r (5.15) 
impUes that \ lm E^+zep+weJ < r. Since 

sinh"^(X) = log (^X + y^X^TTj {MX G M), 

the inequahty (j5.15l) is equivalent to the inequality s < ^ log Ft ^t' ,d{'i'). □ 
We define the perturbed covariance C{n.£,x,y(j)y){zep + we^) as follows. 

L ^ — ' 

ker* 



for any (x, ^, x), (y, 0, y) £ T x {t,i} x [0,^), p,q G {I,-- - ,d} and z,w G C. 
Moreover, if a;,y G [0,/?)^ we write 

Ch(x^2;,y0y)(2ep + we,) = C(x^a;,y0y)(zep + we,). 

The properties of C{zep + weq) are summarized as follows. 

Lemma 5.6. (i) For any (x, ^, a;), (y, 0, y) G T x {t, 4} x [0, ^), p, g G {1, • • • , d}, 

the function {z,w) i— > C(x^a;, y(/)y)(zep + we,) is analytic in the domain 

1 /tt 

(z,w) G |Imz|, |Imw| < -logFt,*.,^ ( - 

(a) The inequality \C{x^x,y(j)y){zep + weq)] < 1 holds for any {x,^,x),{y,(j),y) G 
r X {t, i} X [0, /3), p, g G {I, - ■ ■ ,d} and z,w e C with \ Imz|, | Imui] < 
ilogFt,t,,d(7r/(2/3)). 

(Hi) For any z,w ^ C with \lmz\,\Imw\ < ^ log Ft,t'.d{'^/{'^P)) andp,q^ 

{1, • • • ,d}, the covariance C{zep + weq) satisfies the determinant bound (|5.6p . 

(iv) For any p, 9G{1,-- - ,d},z,w(zC with \ Im z|, | Im w| < ^ log Ft^t' ,d{T^/ P) and 

he2N/p, 

det(C?,(x^x,y02/))(x,j,^),(y,^,j,)grx{t,i}x[o,/?);. = H , f^E^.. -,2 ^ 



32 



Remark 5.7. In (pv)) and in its proof below let us think that each element of 
r X {t,i} X [0, (3)h has been given a number from 1 to N{— 2l3hL'^) and according 
to this numbering the N x N matrix Ch{zep + weq) is defined. However, the claim 
(|iv)) implies that the value of the determinant is independent of how to number the 
elements of T x {t,i} x [0,^)/i. 

Proof of Lemma \5.6\ If | Imz|, | Imwj < ^ \og Ft ^t' ,d{'^ / P), Lemma 15.51 ensures that 
Imi?k+zep+«ieJ < 7r//3. Then, for all k G T* and a G {1, -1} 

[1 _|_ g/3QBk+^ep + „,e, |2 ^ ^g/SaRoBk+^ep+^e, _ j ^ > 

Thus, the denominators in C(x^x, y0y)(zep +weq) do not vanish, which proves 0. 

If |Imz|, |Imw| < i log i^t^t/, d(7r/(2/3)), by LemmaES] | Im£'k+^ep+toeJ < 
7r/(2/3). Thus, the claim dm)) holds true by Proposition 15.31 Moreover, for any 
X e [0, (i) and k e T* 



1 + e^""^^' 



gKRcBk+^ep + u 



- (1 -|_ e2/3Re-E;k+^ep + ™e,)l/2 - 



from which the claim (jn]) follows. 

To show (jiv]) , define an ^.-dependent finite set Mh of Matsubara frequencies by 

Mh := jo; e ^ — ^ I - nh < U! < nh 

Define the unitary matrix Y — 

(y(k0a;,x^a;))(k,0,<^)Gr'x{t,i}xAfh,(x,^,x)erx{t,i}x[o,/3)h by 



lUJX 



Using the assumption that h 6 2N//3, the argument parallel to [71 Appendix C] 
demonstrates that 

{YCuY*WuM^) = ^k.fc^,,j^.,. ^_^,,,/,/^^^^^^^^^^/, . (5.16) 

As in [71 Proposition C.7], the diagonalization (j5.16p deduces (jiv])- □ 

Here we present a calculation, which shows the essence of our analysis to bound 
the correlation functions. Take any n e N and set r„ := ^ logFt_t'_d(/3/(27r)) > 0. 
By periodicity we observe that for any p,q G {1, • • • , d} and z G C with | Imz| < 
ilogFt,«.,d(/?/(27r)) 

e-'-r{^-y^-.)C{^^x,y(lyy){zep) = C(x^a;,y<^y) (ze^^^e^ . (5.17) 
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Moreover, by Cauchy's integral formula we have 



2ti/L 



L f ' d 



2ti Jo 2TTi 7j„_e|=r„ {w - Of 



where dw stands for the contour integral along the contour {w G C | \w — 

(i\ = ^n} oriented counter clock- wise. With this choice of r„, Lemma [5.61 (P allows 
us to repeat this operation n times to obtain 

Lemma 5.8. For any z G C with \ Imz| < i log Ft^t' .d{P/ (Stt)) andp, q ^ {!,■■■ ,d} 
/ e'^~^^~y'^''' — 1 



2t:/L 



C(x^x,y(/)?/)(zep) 



■ C{x£,x, y(py) ^zGp + ^^^65 , , 

w)/iere r„ = -^JogFt^t'Al3/{2'K)). 

Evaluating both sides of (|5.18l) leads to the following bounds. 

Proposition 5.9. For any z G C with \ lmz\ < ^ log Ft ^t' ,d{f3 / {2tt)) , p G {1, ■ • • ,d}, 
the following inequalities hold. 

|C(xea;,y0y)(zep)| <2Ft,t,,rf(^^^ (5.19) 
for all (x, ^, z), (y, 0, y) G F x {t, |} x [0, Especially, 



TT 



|C(xCx,y</)y)(zep)| < 2Ft,t:^, i^ — j (5.20) 

X - y G ({- lL/2\ , - [L/2\ + 1, • ■ • , - lL/2\ +L- IjyK 
Proof. By using Lemma 15.61 (In)) and the inequality 



< 



2n 



log Ft.t'A^/C^f^)) J ' 
(5.21) 
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we have 



I (The right hand side of dEIHl))! < 



2n 



(5.22) 



logFt,t'A{7r/m), 

Stirhng's formula states that for any n E N there exists e{n) e (0, 1) such that 

n" = n!(27rn)-i/2e"-^(")/(i2n) (523) 

(see, e.g, d Eq.6.1.38]). On the convention that 0° = 0! = 1, the formula (15^23)) 
gives the inequality 

n" < n!e" (Vn e N U {0}). (5.24) 
By combining (|5.24p with (|5.22l) we obtain 



1_ / logf,.t'.d(7r/(2/3)) 
n\ V 2e 



27r/L 



|C(x^x,y(/)2/)(^ep)| < 1. (5.25) 



By multiplying both sides of (|5.25p by 2 " and summing over N U {0}, we have 

|C(x^a;,y<^y)(2ep)| <2Ft,t-,rf(^^^ 

for any q G {1, • • • ,d}, which yields (|5.19p . By applying the inequality that |e'^ — 
1| > 2\e\/Tr {ye G [-7r,7r]) to ([5Jg)) . we can derive (jOO)) . □ 

Corollary 5.10. For any z G C iwii/i |Iniz| < ^ log Ft ^t',d{f3/{2TT)) and p G 



l/(2e7rd) 



^e[-/3,/3)hxer 



t,t',d 



(*) 



l/(2e7r(i) 



(5.26) 



Proof. By using (j5.20p and periodicity we have that 
^|C(xex,OC0)(zep)| 



LL/2J+L-1 

|C(xCx,oeO)(zep) 



xer 



vje{i,- - ,£i} 



1=1 



< 2 + 2 ^ F, 
which gives (|5.26p 



t.t'.d 



-2e,d \ / J:'t,t',d [2^ 



l/(2e7rd) 



l/{2e-!Td) 



- 1 



□ 
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Remark 5.11. Remark 1 2. 71 and the inequality (|5.26p imply that 



a;e[-/3,/3)hxer 



as /3 -> +00. 



Remark 5.12. It is the same procedure as the proof of Proposition 15.91 to derive 
the decay bound on the determinant of the covariance. By using the equality that 



'2tt 

^ / / l<j,k<n 



det ( C{xj^jXj,yk(l>kyk) -r^q 



and the determinant bound Proposition 15.31 one can prove that for any (xj , , ) , 

(yj:</'j,%) er X {t,;} x [o,/3) (Vj e {i,--- ,"}) 

|det(C(xj-^jXj,yfc(?!)fcyfc)) l<j,fc<n I 

\ 4e£i 2^5 = 1 27r/i | 



5.3 Proof of the main theorems 

We prepare some lemmas and give the proofs of Theorem 12.41 and Theorem | 
The following lemma fundamentally supports the validity of our argument in this 
subsection. 

Lemma 5.13. For any Rr > 0, Ri E (0, | log Ft.f'_c((7r//?)) and p E {!, - ■ ■ ,d} there 
exists Q > such that the function 

{z,T]) I — > (C/i(zep),?7) 

is analytic in the domain 

{(z,7]) e I |Rez| < Rr, \lmz\ < Ri, \r]\ < Q}. 

Proof. Note that by Lemma [5^ (HI) for any Rr > and Ri G (0, ^ log Ft ^t' Att/ (3)) 
the function z C;i(x^a;, y(j}y){zep) is bounded in the compact set G C | | Re z\ < 
Rr, |Imz| < Ri}. Thus, by Lemma [ 

lim sup 

V-^O zee 
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Hence, we can take a small Q > such that the denominator of Sj^ri^ y-m i,™ 
(Ch(zep), rf) does not vanish for any rj € C with \r]\ < Q and z G C with | Re z| < Rr, 
lmz\ < Ri. This proves the analyticity of 5^^,,, y,,, ~a 4™ (C'/i(zep), 77) in the same 
domain. □ 

For p G {1, ■ ■ ■ ,d} we define the matrix Up — {Up{{'x.^x)j , {^£,x)k))i<j,k<N by 

The following equality is based on the C/(l)-invariance property of the Grassmann 
integral J ■dipxdipx- 

Lemma 5.14. For any p e {l,---,d}, f{:ijjx,^x) £ C[V'(x^^)^. , V'(x4:c),- I j e 
{1, ■ • • , A^}] and z £ C with |Imz| < ^ log Ft_t/^d(7r/(2/3)), the following equality 
holds. 



fitPx,^X )dlJ.C^ {ze^)ii^X,^x)- 

Proof. Observe the invariance that 

f{m^JxYAl^P^xy)d^xdJjx = J fi^x.Md^xd^x, 
which implies that 

fiiU^i^xY, (5.27) 



X- 



(5.28) 



Moreover, by using the equality (j5.17p we have 

{Upipl(,Ch{zepy^Up^x) = {i^xd^pCh{zep)U^y^i^x) 

= {i;'x.Ch{{z + 2n/L)ep)-^'x)- 
By substituting (|5.28p into (|5.27p . we obtain the desired equality. □ 
Combining Lemma 15.131 with Lemma l5 . 141 shows 
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Lemma 5.15. For any n G N, set fn ^ log Ft^t' ^i"^/ (2/3))- Let p G {1, • • • ,d}. 
There exists Qn > such that for any 77 G C with \r]\ < Qn, 



2tt/L 



where §^^._g.^^^ dzj stands for the contour integral along the contour {zj G C | \zj — 
^il = ''n} oriented counter clock-wise. 

Proof Take any R, G (i logFt,t,,rf(^/(2/3)), i log i^t,t^d(7r//3)) and i?, > nr„ + 

By Lemma 15.131 we can take a small Qn > so that (z, r/) 1— S'j^™ y^, 

(Ch(zep), 77) is analytic in the domain {(z, 77) G | |Rez| < Rr, |Imz| < Ri, \r]\ < 

Qn}- 

Fix any 77 G C with jryj < Qn- By noting this domain of analyticity, Lemma 
15.141 the fundamental theorem of calculus and Cauchy's integral theorem verify the 
following equalities. 



2-k/L 

L /-^"/^ ,„ 1 /■ , 1 



(5.30) 

Then, we multiply both sides of ((OD)) by (e'T(E7=iij-E"Liy..ep> _ i)/(27r/L) 
and do the same calculation for the integrand Sj^rh y™ i,™ (C^ (zie^) ,ry). Re- 
peating this procedure 77 times results in (j5.29l) . □ 

Lemma 5.16. For any 777,71 G NU {0} let bn,m denote the coefficient of rf^ in the 
Taylor series of the function 



2tt/L 
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around G C. The following bounds hold. For any n G N U {0}, m G N, 



|6„.o| <2"n" logi^, 



1 

m 



t,t'.d 



t,t' ,d 



Ft,t'd 



2(3 
n 



to16" 

l/(2e7rd) 



l/(2e7rd) 



1=1 



J 



where 0*^ = 1. 



Proof. We apply Proposition 14.11 to bound the integrand 
'5'x™,y™,H™,4.-(C'»C"=i^jep),7]) in (|5.29p. Since jlmXl^Li^il < 
i logFt_t/_d(7r/(2/3)), Lemma (Iiv)) imply that the covariance C;i(^"^^ ^j^p) 
satisfies the assumptions of Proposition 14.11 with B = 4. Moreover, Corollary 15. 10[ 
the periodicity and the translation invariance of Ch ensure that 



V <A(3 



Ft 



l/(2e7rd) 



l/(2e7r£i) 



Therefore, if bm denotes the coefficient of ry™ (Vm G N U {0}) in the Taylor series 
of the function 77 1— > 5"^^^, Srsi <|m (C'/i(X]j=i ^3^p)tV)i obtain the inequalities 
that |6o| < 4™, 



/ 



\b,n\ < 



ml6" 



/3 



Ftf, 



l/(2e7rd) 



+ 1 



Ftf, 



l/(2e7rd) 



E'16'l|f/L,d|L. 



1=1 



(Vm G N). 



(5.31) 

In the right hand side of (|5.29l) we expand Sj^^ $™(C?i(X]J=i ^j^p)^v) 

into the power series of 77. Here note that by Fubini-Tonelli's theorem we may 



exchange the integral H^Li /q dOj /| 

m=0 lli=l \2iT Jo 

thus the equality 



/L 



d0 



dzi 



27r/L 



if 1 77 1 is sufhciently small, and 



2^ Jo 



d6. 



1 



dzi 



holds for all m G N U {0}. Combining the inequality (I5.2ip with the bounds (|5.3ip 
yields the desired bound on |6„.m|. □ 
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Corollary 5.17. Assume (j2.8p with R G (0, 1) and use the same notation as in 
Lemma \5.16l For any n G NU{0} and p G {1, • • • ,d} there exists Nq G N such that 



2tt/L 



(5.32) 



( X) j = 1 S j = 1 yi ' 



2tt/L 
<2"n'^ [log Ft^t,^d 



(4™ -ml6™log(l -i?)) , 



(5.33) 



for all he2N/P with h > 2Na/l3. 

Proof. Since Tre~'^(^"+''^)/Tre~'^^" > for all 7y G M, the uniform convergence 
property (j3.13l) ensures that there exists iVo G N such that 



> 
(5.34) 

for all 7^ G M with \ri\ < 1 and h G 2N//3 with h > 2No//3. Let us fix such a large 
h. Since the Grassmann Gaussian integral in (|5.34l) is a polynomial of 77, we can 
take a domain Oh C C such that [—1,1] C Oh and the inequality (|5.34l) holds for 
all 77 G Oh- This proves that the function 



?7 



27r/L 



is analytic in the domain Oh- 

On the other hand, Lemma [5.16l and the inequality ||?7L,il|L,; < \\Ui\\i imply that 
if 7] G C satisfies 



/ 



\v\ < 



Ft,f,d ( 



l/{2eiTd) 



F, 



l/{2eTvd) 



E^16'||C/H| 



- 1 / i=i 



i?, (5.35) 



E IV™lh'"l<2"n" logFt,*,,, ^ 



2"n" logF,,,.. — 



2/3 



00 

116™ y — i?" 

' ^ m 



2"n" logFt^i,,^ 



(4™ - ml6" log (1 - i?)) < cx). 

(5.36) 
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Since by the assumption (j2.8p the right hand side of (|5.35p is larger than 1, the 
identity theorem for analytic functions proves the equality 



1 



2tt/L 



for ah ?7 e [-1, 1]. The bound (jOSI for = 1 gives 
As the last lemma let us claim 



□ 



Lemma 5.18. Assume (j2.8p . The thermodynamic limit (|2.9p exists and takes a 
finite value. 

For continuity of our argument in this subsection we present the proof of Lemma 
15.181 in Appendix [B] We now proceed to 

Proof of Theorem \2.4\ By substituting (j5.24p into (j5.33|) and dividing both sides by 
4"e"n! • (log Ft,t/,d (7r/(2/3)))"" we have 



1 /l0gi^M',<i(§) 



4e 



2tt/L 



•Sjj^m^ym em I,™ {Ch, 1) 



< 2-" (4™ - ml6" log(l - R)) 



(5.37) 



for aU 71 e N U {0} and p e {!,■■■ ,d}. Summing (|5.37p over n e N U {0} yields 



S 



,^ACh,1) <2(4™-ml6™log(l-i?)) 



1^ sr-'d 

^ \ 4ed Z^p = l 



27r/i 



(5.38) 



for all /i S 2N//3 with h > 2No//3. To derive the same bound on Sym |,m i. 
(C/i,l) from (|5.38p is immediate. Then, Proposition 13.91 ensures that 



xiji ^x^Cri, ^ymipm ^yiipi ■ ^yi0i 
< (4^+1 _ m42"+i 



1^ v^d 

^ \ 4ed Z^p=l 



log(l - R)) 



Finally, by Lemma [5. 181 we can take the limit L — > +oo in the both sides to obtain 
the claimed inequality. □ 
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Proof of Theorem \2.5\ Since the construction of Theorem 12.51 is parallel to that 
of Theorem 12.41 we only outline the proof. The main difference is that we use 
Proposition 14. 21 in place of Proposition 14. II 

(Step 1) For any m,ri G N U {0} let Cn,m denote the coefficient of 77'" in the 
Taylor expansion of 



gj^{xi+x2-yi-y2,ep) „ Y 
2n/L 



(ii,i2),(yi,y2),(t.i).(t.i) 



around € C. By repeating the same argument as in Lemma lS. 161 using Proposition 
4.21 in stead of Proposition 14.11 we obtain 



<2"7i" logi^^t,*,,. 



2/3 



64 



3m + 4 
3m + 4 V m 



16/3 



\ 



l/{2eTvd) 



I 



^ Vi^M',d(^) -1^ 
(Step 2) Recall the fact that the radius of convergence of the power series 
^' 4 / 3m + 4 



E 

m=0 



3m + 4 



is 4/27 and 



E 

T?l = 



3m + 4 



3m + 4 



81 
16 



(see [21 Lemma 4.9]). 

Therefore, the argument involving the identity theorem parallel to the proof of 
Corollarv 15 . 1 71 shows that for all C/ e M satisfying 



\V\ < (108/3)-i 



t.t'.d 



(*) 



l/(2e7rd) 



F 



t,t' ,d 



2/3 



l/{2e-iTd) 



and for sufficiently large h G 2N/(3 the equality 



Cgj^(xi+X2-yi-y2.ep> _ ^\ 

and the inequality 

gi^(xi+X2-yi-y2,ep) _ ^ 



S'(xi,x2),(yi.y2).(t:i):(t,i)(^'" 1) - E 



27r/L 



<2V [log Ft,t',d 



2/3 



'(ii.X2),(yi.y2),(t,i),(t,i) 

— n 

•81 



(5.39) 
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hold. By using the fact that S'(xj.x2).(yi,y2),(t:^):(t4) (^'" ^) ^ rational function of 
U, we can prove (jOO)) for f7 6 M with ((^TT|) as well. 

(Step 3) By repeating the same calculation as in the proof of Theorem 12.41 or 
Proposition 15.91 we reach the bound 

1 I ^i2,r(ici+x2-yi-y2.ep>/-L_^ | ( K ACW 



< 324 • Ft.t'.d 



2/3 



for [/ e M satisfying (pUj) from ((0^)) . 
(Step 4) By writing 

U Ct^ii^xi^xt 
xer 

= C^o((xi,X2),(Cl, 6), (</■!, '/'2))V^;i5iC26^''202V'xi0i 

xi,x2er5i,52,0i,02e{t4} 

with C/o((xi,X2), (Ci^6), (01,02)) := C^'^xi,x2'^5i,T%4'^«i,0i%,02 7 one can tr anslat e 
the proof of Lemma lS.lSl to confirm the existence of the thermodynamic limit (j2.12D . 
Then, by sending L +oo in (j5.40l) we obtain the claimed bound for U satisfying 
(|2ll|l . □ 

A Anti-symmetrization of the coefficients 

For simplicity we use the notations X', e T', S', e {t,l}' {1 G N) defined in 
the subsection 13. II Moreover, let Xk, Y| denote 

((xi,a),(x2,6),-- - ,(x/,6)), ((yi,0i),(y2,02),-- - ,(yi>0) e (r x {t,;})', 

respectively. Let Si be the set of permutations over I elements. For tt, r G 5*2 we 
write 

'"'(^k) ((X,r(l),^7r(l)), • • ■ , (X7r(/),^7r(0))' 
^(YI) ■■= ((yT(l),0T(l)), • • • , (yr(/)>r(;)))- 

Lemma A.l. Assume that I e {!,■■■ ,21"^}. For any function g : x {t, lY C 
there uniquely exists a function / : (F x , 4-})^' C such that 

finixL),T{Yl)) = sgn(^) sgn(T)/(X|, 4) (A.l) 

/or aZ/ TT,T e Si, XL,Yl e {T X {t, i})' and 

I] Y 5(^',S',$')V';,j^ •••i^i^j^Vx,^, •••V'xi^i 
x'er' H'.*ie{t.4.}' 

(A.2) 

^|,>1,6(rx{t,;})' 
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Moreover, the following inequality holds. 



maxl max V V f{XL,Yl)l 



(rxlt.lTV (yfc.*fc)erx{T,4.} J 



max 

(yi,0i)erx{t,i} 

" VfcG{2,---,i} 

< max < max max lofX', S', $')|, 



(A.3) 



max max 
je{i.--- ,1} (xi.<hj)erxH,i} 

vfce{i,... ^' 

Proo/. By setting u{Xl,Yl) := S', we see that 

E •••'0;,{>yi0i •••^y,0, 

-Y^Xieirxft.i})' 

-Y^X^eirxft.i})' ^^,^65, ■'' 
>xi{i ■•■V'x,6^yi'>i ■■•^yi</'i- 
If we define / : (r x {t, ^ C by 

f(^ri vi\ Sr sgn(7r) sgn(T) , , 

f satisfies (|XT|) -(|X2 |) and 

max y y \f{xl,Yl)\ 

(xi,6)erx{t,i} ^ ^ \J\ 4-^1 

< max max y y \u{X~,YI)\ , 

' ' (-fe.Sfe)erx{t,4} yig(rx{t,4.})' ^ ' 



max max " " i^i v- 



vfce{i,.-. ,i}\{j} ^' 
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The inequality 

max y y \f{xi,Yi)\ 

^^^^^'■■■■'^ (A.5) 
< max max V V S', 

je{i,-,/}(x,,0,)erx{t,i} , ^ t^, 

vfce{i,... 

can be confirmed similarly. The inequalities (|A.4I) - (IA.5I) yield (jA.3|) . 

To prove the uniqueness of such /, we number each element of F x {t, 1} so that 
r X {t,|} = {(x^)j I j = 1, • • • , 2L'^}. Suppose that 

J2 /(^L4)V';iCi •••^;,?,^yi^i ••■'^y^^i =0- 



The anti-symmetricity (jA.ip implies that 



l<31<-<3i<21-'' 
l<fel<-<fej<2i'^ 



(A.6) 



•V'(x5),, •••V'(x«),,V'(xS)., •••'/'(xc),, =0. 

Since 

{■^(■xOji ■ • • V'(*xC),^V'(x5)fc^ ■ ■ • V'(x5)fcJl<ji<...<j,<2L'i, l<fci<---<fc,<2L'l 

are linearly independent in the complex vector space of linear operators on the 
Fermionic Fock space on F x {t, 4-}, the equality (IA.6I) deduces that 

fii^On > • • • , i^On : (xOfei , • • • , (xOfc, ) = (A.7) 

for all jm, km e {!,■■■ , 2i'^} (m e {1, • • • , /}) with ji < ■ ■ ■ < ji, h < ■ ■ ■ < h. 
Again by ((XT|) the equality (|X7)) ensures that /(A:4,rj) = for all Xi^Y^ e 
(F X {t,4-})', which concludes the uniqueness of / satisfying (jA.ip and (|A.2p . □ 

B Proof of Lemma 15.181 

By Corollary 15 . 1 71 with the same notation there exists Aq G N such that 

oo 

%-,y-,H-,4.-(^'- 1) = E ^o,m (B.l) 

m=0 

for all h G 2N/^ with h > 2A^o//3- By LemmatSUl] |6o.o| < 4" and 

rr? 16™ 

|6o ml < ^-^i?" (Vm e N). (B.2) 

TO 
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Assume that the hmits 

Um 6o.m, hm hm 6o.m (B.3) 

exist for any m G N U {0}. Since the right hand side of the inequahty (|B.2p is 
summable over to G N, we can apply the dominated convergence theorem for Z^(NU 
{0}) to verify that the equalities 



lim 6o.m = hm 6o,m, 

he2N//3 m=0 m=0 fee2N/3 

oo oo 

lim lim y &o m = > hm lim 6n . 

ieN )ie2N/3 m=0 m=0 -C-gn )ie2N/3 



(B.4) 



hold in C. Combining (jB.4[) with (jB.ll) proves the existence of 

lim lirn 5^^ |,™(Cft, 1) 



for ah (X™, y™, S", $") G (Z'^)2™ X {t, 4-}^'". Proposition [SJ then concludes that 
the correlation function converges to a finite value as L — +cxd. Hence it suffices 
to prove the existence of the limits (IB.3p . 

By (021) for Gh = Ck, 6o,o = det(C(xj4o, yfe'^fcO))i<i,fc<m- From it is 

clear that limi_>+oo,LeN &o,o exists. 

Let us derive an expression of 6o,m (jn > 1) from the tree formula (|4.5p . We 
consider that the root of any tree T G T({0, 1, • • ■ , to}) is the vertex 0. For any 
j G {0, • ■ • ,to} let disrO) (G {0, • • • ,to}) denote the distance between the root 
and the vertex j along the unique path of T connecting with j. We write any 
line of T as {p,q) {p,q G {0, • • • ,™}) with disr(q) — disT(p) + 1- Set max(r) :— 
maXjg{i_..._„}disT(j) and Dj{T) := {{p,q) G T \ Aist{p) + 1 = disT(q) = j] for 
j G {!,••• ,max(r)}. 

By letting J|^^ ^■^^j,{Ap^q + Aq^p) act on the monomial (|4.7p first and then ap- 
plying e'^^-^-^'^'^''^'^^^ to the rest of the monomials in (j4.5p . we have 



bo,. 



E E E ^^/:'-'"""°''(-.^oj 

JO,g-lfco,g-lao,q€{-l,l} 

max(T) / /p 

• n n E E E 

r=2 (■p,q)£Dr(T) \jp,5 = lfcp,, = lap,,e{-l,l} 
' f{T^ {jp,q7 kp,q> ^P,q} {p,q)eT j ^h), 

(B.5) 
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where for (0,(j) S Di{T) 

(^Jo.,:feo,,,i(^^^ := C'/i(x^g,feo,,S9:yjo.,'^jo,,So), 
for ip,q)eDr{T) (r >2) 

(x,y) := C„(xep,,,„sp,y0,,fc^,,s,), 

and the term /(T, fcp^g, o,p,q}{p,q)£T^ Ch) satisfies that if there exist (w, q), (u, g') 
e T such that q ^ q' and {jv,q,av,q) = ijv,q' jO-v.q') or if there exist (p,v), {v,q) G 
T such that {kp^y,ap^y) = {jv,q, -av,q), then fiT,{jp^q,kp^q,ap^q}(^p^q)(zT,Ch) = 0. 
All the C;i-dependent terms contained in f{T,{jp^q,kp^q,ap^q}(^p^q'j(zT,Ch) are the 
remains left after e^'^''^^-'"''^'^)) acting on Grassmann monomials. 
Moreover, by setting 

f'{T, {jp^q, kp^q, ap,g}(p^q)gT, Ch, {ULdk}k=l) 

m 

■■= f{T, {jp,q, kp,q, ap,J(p.,)eT, Ch) J] {X^ , , $1'= , 

k=l 



the expression (jB.sp is rewritten as follows. 
bo,m = 

I 



soe[0,/3)fc k=l 



\ 



??7 



E 

TeT({o,i,--- 



(0,g)GDi(T) \ Jo,q-lfco,g-lao,qG{-l,l} 

max{T) / Ip Iq 

n n E E E 

r=2 {p^q)(^D^(T) \ Jp,5 = lA:p,, = lap,,e{-l,l} 



• C'/i"'" (%Jp,, 'Xg.fcp.J J /'(^; {jp,q,kp^q,ap^q}(^p^q)£T,Ch, {t^L.ij, jfe^l ) ■ 

(B.6) 

Note that by Proposition 15.31 the bound of the form (|4.8p for ^? = 4 is valid. By 
using this bound, the equality (14.61) and the inequality that 



47 



we can find a non-negative function g(T, {jp^q, kp^q, ap^q}(p^q)^Ti {lk}]^=i), wliicli is 
independent of L, h and the variables , S'j!" , ^J.*")}^]^, such that 

\f'{T,{jp,q,kp,q,ap^q}(p^q)(zT,Ch,{UL.lkyk=l)\ 
< g{T, {jp,q, kp^q, ap,q}(p,q)eT: {^felfcLl)- 

By using the periodicity and the translation invariance of Ch and Ul.i and by 
arguing recursively with respect to r running from niax(T) to 1, we see that 



n 

fc=i 



E 



(rh 
Jo,g-lA;o,q-lao,qG{-l,l} 
max(T) / Ip Iq 

n n (uL,Mi^,^[^,K^y'^>-^ E E E 

r=2 {p,q)£D^{T) \ Jp,9 = l fcp,9 = l ap,<,e{-l,l} 

' ^h""' ""' ^"(^PJp.g' ^9>fep,g)^ f'{T, {jp,q, fcp,q, ap,q}(p,g)GTj C*/! j {f^i./fclfcLl) 
ml \ / ni 

-H E n EE E ^^r"'-"-""'K.oj 



k=l 



■ UL,i^{{^qs,- ■ ■ ,X5,fc„_,_i,X5,fc„_^+i, • • • ,x,j,, 0), j 

max(T) / /p 

• n n E E E ^^^--'-'"-'(o,x,,,^j 

i-=2 (p,<j)G-Dr(T) \jp,5 = lfcp,, = lap,,G{-l,l} 

• f/L,i,((x,,i> • • • ,Xg,fcp,,-i,x,,fc^ ,+1, • • • ,x<,,,,,0),S^',$^«)i'«s=^ j 

• f'{T, {jp^q, kp^q, ap^q}(^p^q^(zT, Ch, {UL,h}T=l), 

where 

(x) := C;,(o4„so,x0,,fc„„s,), 
Cf-'=°-\x) C;,(x^,,fc„„s,,O0,„„so). 

Remark that the dependency of /'(T, {jp^^, kp^q, ap,q]{p,q)(^T, Ch-, {Ulj^]'^^i) on the 
variables {X^^ , • • • , ) is changed in every step of the derivation of (|B.8|) . However 
we used the same notation for simplicity. 



(B.8) 
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By substituting (|B.8|) into (|B.6|) we obtain 



HE E 



m n 



fc=l 



(B.9) 



L.h 



where 



<'.'i<.)i-..«^!''-i''*'i'''----<^''"-="'*"'' 



fc=0 \ Sfce[0,fi)fe/ A;=l 

EE E 

> - , io.5 = lfco,, = lao.,e{-14} 



<'"6(({-LV2J,-,-Li/2j+i-i})<^)''= 



E 



TeT({0,l,--- ,ni}) 



jo,g ifco.g ,ao,5 



■ n f E E E ^^^''■"''"'"^''(o.x,.,,,) 

(p,q)eT\Di(T) \ip,,=lfep,, = l ap., £{-1,1} 

• ?7/,((x,,i,-- - ,x,,fe^^_i,x,^fe^^+i,-- - ,x,,,,,0),S^',$^')i'<j^2j 
By the properties of Ch and the equahty (|2.5p for ^ = 1 , we observe that the Hmits 



exist. 



lim i? r, im , hm hm B ^m. 

he2N//3 Z.eN /ie2N//3 



By using the inequalities (|5.20l) and (jB.7[) we have 



(B.IO) 
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with 



Clp Iq 
E E E 
_ , , , ,, ip.5=ifep,<j=iap,,e{-i,i} 



5(2^5 flp,g}(p,g)GTj {^fe}fcLl), 



where /q i^- Note that i?^ ,r;j.i,m_^ is independent of L and /i. Moreover, by 
calculating in the same way as in the proof of Corollary 15.101 we obtain the bound 
that for any m G N 



m I n 

HE E 



{XI'' ,3]!' ,<i>i'')e{i.'')'kx{tA} 



( 



< 



n 

k=l 



2(3 



t,t'd 



Clp Iq 
EE E )5(?',{jp,g,fcp,g,ap,g}(p,5)eT,{'fc}r=l) < 

(B.U) 



(ft)'"' 




1 / (^=1 



TeT({0,lr-- 



By (IB.gp - dB.lip we can apply the dominant convergence theorem to ensure that 
the limits (IB.Sp exist and satisfy 



m I 71 



Jim^ bo,m = n E E 



lim B 



L.h 



/ 



lim lim &o,m = TT E E 

fc=i y;,=i (x^'=,H^^<i>'^'=)e(Z'')'fc x{t,i}"'fc / 



hm lim B r, i™ , 

iSN (,g2N//3 



which completes the proof. 
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